Graduate Advisor
April 27th, 2018
object due to ~ symbol => impute ~ as NaN and parse data type to float[=] check duplicates: explain pandas table, duplicated galaxies with Dist > 30 Mpc => drop
extract local galaxies within 30 Mpc and 3D view
spatial distribution and edge galaxies
check interesting features
Description
1- 7 I7 --- PGC [2,4715229]? Identifier from HYPERLEDA
(empty for globular clusters)
9- 36 A28 --- Name Common name of galaxy or globular
38- 46 F9.5 h RAhour Right ascension (J2000, decimal hours)
48- 55 F8.4 deg DEdeg Declination (J2000)
57- 60 F4.1 --- TT [-9,10]? Morphological type code (1)
62- 66 F5.2 mag Bmag ? Apparent blue magnitude
68- 74 F7.3 arcmin a ? Major diameter (arcmins)
76- 82 F7.3 arcmin e_a ? Error in major diameter (arcmins)
84- 90 F7.3 arcmin b ? Minor diameter (arcmins)
92- 98 F7.3 arcmin e_b ? Error in minor diameter (arcmins)
100-104 F5.3 --- b/a [0,1]? Ratio of minor to major diameters
106-110 F5.3 --- e_b/a ? Error on ratio of minor to major diameters
112-116 F5.1 deg PA [0,180]? Position angle of galaxy
(degrees from north through east)
118-123 F6.2 mag BMAG ? Absolute blue magnitude
125-131 F7.2 Mpc Dist ? Distance (Mpc)
133-138 F6.2 Mpc e_Dist ? error on Distance (Mpc)
140-143 F4.2 mag e_Bmag ? error on Apparent blue magnitude
145-148 F4.2 mag e_BMAG ? error on Absolute blue magnitude
'object' (count 19): - Name - Type - App_Mag_B - err_App_Mag_B - Abs_Mag_B - err_Abs_Mag_B - App_Mag_I - err_App_Mag_I - Abs_Mag_I - err_Abs_Mag_I - App_Mag_K - err_App_Mag_K - Maj_Diam(a) - err_Maj_Diam - Min_Diam(b) - err_Min_Diam - b/a - err_b/a - PA 'float64' (count 4): - RA - Dec - Dist - err_Dist 'int64' (count 1): - PGC
| Name | Dist | err_Dist | RA | Dec | Type | App_Mag_B | |
|---|---|---|---|---|---|---|---|
| 63638 | PGC138606 | 34.347 | 5.152 | 3.18005 | 61.105202 | 5.0 | 16.32 |
| 68443 | PGC166081 | 72.361 | 15.919 | 4.09247 | 71.469597 | 7.0 | 17.47 |
| 68880 | PGC166081 | 24.847 | 3.727 | 4.09246 | 71.469470 | NaN | 17.39 |
| 130474 | PGC138606 | 35.569 | 7.825 | 3.17891 | 61.113250 | 1.2 | 15.51 |
| 178917 | 6dFJ1705055-200214 | 112.653 | 24.784 | 17.08485 | -20.037250 | NaN | NaN |
| 178919 | 6dFJ1704153-203840 | 117.694 | 25.893 | 17.07090 | -20.644350 | NaN | NaN |
| 179298 | 6dFJ1705055-200214 | 112.653 | 24.784 | 17.08485 | -20.037250 | NaN | NaN |
| 179620 | 6dFJ1704153-203840 | 116.194 | 25.563 | 17.07090 | -20.644350 | NaN | NaN |
[x] check duplicates
[=] extract local galaxies within 30 Mpc and 3D view: 8,946 galaxies out of 183,327, present in 3D plot
spatial distribution and edge galaxies: present in distplot
check interesting features
[x] check duplicates
[x] extract local galaxies within 30 Mpc and 3D view
[=] spatial distribution and edge galaxies: present in distplot
be cautious with bin in histgrams[-4, -3]: E/S0 => GCpG.iMType==6[0, 11]: S/Irr => GCpG.iMType==1[-6, -4): E => GCpG.iMType==2(-3, 0): S0 => GCpG.iMType==3NaN: Unknown, 0[x] check duplicates
[x] extract local galaxies within 30 Mpc and 3D view
[x] spatial distribution and edge galaxies
[=] check interesting features: missing data
| Name | RA | Dec | Dist | Type | Abs_Mag_B | Abs_Mag_I | |
|---|---|---|---|---|---|---|---|
| count | 8946 | 8946.000000 | 8946.000000 | 8946.000000 | 5979.000000 | 7877.000000 | 5606.000000 |
| unique | 8946 | NaN | NaN | NaN | NaN | NaN | NaN |
| min | NaN | 0.005940 | -89.334590 | 0.010000 | -5.000000 | -22.590000 | -25.520000 |
| max | NaN | 23.996680 | 86.131800 | 29.986000 | 10.000000 | -4.970000 | -1.680000 |
| mean | NaN | 10.973250 | 3.455383 | 19.011495 | 4.524034 | -16.210734 | -17.158776 |
| std | NaN | 5.277991 | 34.889882 | 6.616150 | 4.707900 | 2.451796 | 2.627928 |
Description
0 1- 14 A14 --- Name Galaxy identifier
16- 27 A12 --- OName Other galaxy identifier
29- 38 F10.7 h RAhour ? Right Ascension in decimal hours (J2000)
40- 50 F11.7 deg DEdeg ? Declination in decimal degrees (J2000)
4 52- 56 A5 --- MType Morphology type
5 58- 63 F6.2 Mpc Dist ? NED distance; see Section 2
65- 69 F5.2 Mpc e_Dist ? Uncertainty in Dist
71- 83 A13 --- n_Dist Method used to determine Dist
85- 89 F5.3 mag Av NED foreground V band extinction
9 91- 97 F7.3 mag VMag [-25/-11]? Absolute visual magnitude (1)
99-104 F6.3 mag e_VMag [0.2/1.5]? Uncertainty in VMag
11 106-115 F10.3 mag KMag [-28/-14]? Absolute K band magnitude (2)
117-125 F9.3 mag e_KMag [0.1/0.3]? Uncertainty in KMag
127 A1 --- l_Ngc [>] Limit flag on Ngc
14 128-134 F7.1 --- Ngc [0/32500] Number of Globular Clusters (NGC)
136-142 F7.1 --- e_Ngc [0/20000] Uncertainty in Ngc
144-152 A9 --- r_Ngc Reference(s) for Ngc (see refs.dat file)
17 154-159 F6.2 km/s sigma [10/414]? Stellar velocity dispersion (3)
161-165 F5.2 km/s e_sigma ? Uncertainty in sigma
19 167-171 F5.2 kpc Reff [0.1/55]? Effective radius (4)
173-177 F5.2 kpc e_Reff ? Uncertainty in Reff
21 179-184 F6.3 [Msun] logMd [7.6/12.8]? Log of dynamical mass
186-191 F6.3 [Msun] e_logMd ? Uncertainty in logMd
23 193-197 F5.2 [Msun] logMGC ? Log of total globular cluster mass (5)
199-203 F5.2 [Msun] e_logMGC ? Uncertainty in logMGC
25 205-209 F5.2 [Msun] logMB ? Log of central supermassive black hole mass(6)
211-215 F5.2 [Msun] E_logMB ? Upper uncertainty in logMB
217-221 F5.2 [Msun] e_logMB ? Lower uncertainty in logMB
MType: str => to numeric, & parse in Harris manner, present in barplot[=] check duplicates and NaNs: explain pandas table, 418 galaxy enteries
[VMag, KMag, Dist], MilkyWay => drop[VMag, KMag], A1689-BCG, 770 Mpc => droplogMGC, 3 in Ngc => data is dirty, check in 4-A-bcheck interested features
Our catalog contains 248 ellipticals, 93 S0's, and 81 spirals or irregulars, (Harris et. al. 2013).
| Name | RAhour | DEdeg | Dist | e_Dist | VMag | e_VMag | KMag | logMGC | Ngc | r_Ngc | sigma | Reff | logMd | logMB | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | MilkyWay | NaN | NaN | NaN | NaN | -21.30 | 0.3 | NaN | 7.66 | 160.0 | 44 | 105.0 | 0.70 | 9.856 | 6.61 |
| 356 | A1689-BCG | 13.191528 | -1.338056 | 790.00 | 40.00 | NaN | NaN | NaN | 10.10 | 30000.0 | 66b | NaN | NaN | NaN | NaN |
| 227 | NGC4417 | 12.447387 | 9.584243 | 16.03 | 0.53 | -20.00 | 0.2 | -22.860 | 7.24 | 72.0 | 77 | 135.2 | 1.22 | 10.318 | NaN |
| 228 | NGC4417 | 12.447387 | 9.584243 | 16.03 | 0.53 | -20.00 | 0.2 | -22.860 | 7.24 | 72.0 | 77 | 135.2 | 1.22 | 10.318 | NaN |
| 272 | VCC-1386 | 12.530939 | 12.656659 | 16.00 | 2.00 | -17.34 | 0.2 | NaN | 6.67 | 26.6 | 67 | 20.5 | 1.80 | 8.848 | NaN |
| 273 | VCC-1386 | 12.530939 | 12.656659 | 16.10 | 1.13 | -16.90 | 0.2 | NaN | 6.64 | 26.0 | 24 | NaN | NaN | NaN | NaN |
| 138 | NGC2915 | 9.436528 | -76.626389 | 4.06 | 0.20 | -16.12 | 0.2 | -18.298 | NaN | 3.0 | 66a | NaN | 0.31 | NaN | NaN |
[x] check duplicates and NaNs
[=] check interested features:
[=] statistics, present in 2 pandas tables
var(Ngc(6:0)) significantly larger => add logNgc(6:-inf)Ngc, e_Ngc(0:3), logMGC(1:NaN), e_logMGC(2:0): Ngc derived from logMGC => focus on Ngc/logNgc(412/418)Ngc == 0: [7,51,87,88,271,304] => Ngc = 0.5 e_Ngc, avoid -inf temporatelye_Ngc == 0: [1,28,138] => e_Ngc = 0.5 Ngc, avoid -inf temporatelylogMGC, e_logMGC == NaN: [138] => linear interpolate NaN logMGC from logNgc(412))e_logMGC ==0: [1,28] => leave alonevisualization
| Name | VMag | KMag | Ngc | e_Ngc | r_Ngc | sigma | Reff | logMd | logMGC | Dist | logMB | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| count | 418 | 418.000000 | 344.000000 | 418.000000 | 418.000000 | 418 | 271.000000 | 340.000000 | 253.000000 | 417.000000 | 418.000000 | 64.000000 |
| unique | 418 | NaN | NaN | NaN | NaN | 106 | NaN | NaN | NaN | NaN | NaN | NaN |
| min | NaN | -24.190000 | -27.080000 | 0.000000 | 0.000000 | NaN | 10.400000 | 0.130000 | 7.693000 | 4.880000 | 0.030000 | 0.000000 |
| max | NaN | -11.170000 | -14.510000 | 32500.000000 | 13000.000000 | NaN | 414.000000 | 55.000000 | 12.726000 | 10.120000 | 375.000000 | 10.320000 |
| mean | NaN | -19.127943 | -22.583840 | 1512.699761 | 351.961962 | NaN | 168.750923 | 4.287588 | 10.727692 | 7.369305 | 34.758971 | 5.933750 |
| std | NaN | 2.799122 | 2.659367 | 3957.436127 | 1075.182192 | NaN | 93.365634 | 6.092928 | 0.969581 | 1.229806 | 46.445852 | 3.953719 |
| Name | Ngc | e_Ngc | logMGC | e_logMGC | |
|---|---|---|---|---|---|
| 7 | NGC221 | 0.0 | 1.0 | 6.26 | 0.10 |
| 51 | FCC-0059 | 0.0 | 4.0 | 6.21 | 0.19 |
| 87 | FCC-242 | 0.0 | 3.3 | 5.51 | 0.30 |
| 88 | FCC-246 | 0.0 | 3.4 | 5.51 | 0.30 |
| 138 | NGC2915 | 3.0 | 0.0 | NaN | NaN |
| 271 | VCC-1363 | 0.0 | 3.0 | 6.19 | 0.24 |
| 304 | VCC-1714 | 0.0 | 4.5 | 6.19 | 0.16 |
[x] check duplicates and NaNs
[=] check interested features
[x] statistics, present in 2 pandas tables
[=] visualization
logNgc: present KDE in distplot NaN data: present in msnoplot[=] explore correlations for data type in float: present in pairplots with heatmap
'Dist', 'e_Dist', 'Av', 'VMag', 'e_VMag', 'KMag', 'e_KMag', 'Ngc', 'e_Ngc', 'sigma', 'e_sigma', 'Reff', 'e_Reff', 'logMd', 'e_logMd', 'logMGC', 'e_logMGC','logMB', 'E_logMB', 'e_logMB', 'logNgc'xxx & e_xxx: intuitive, not interestingAv & Dist: Dist derived from NED Av => focus on Dist (more intuitive)VMag & KMag(74:NaN): intuitive, LR in 2-ADist & Ngc(6:0) & Reff(78:NaN)VMag & KMag(74:NaN) & sigma(147:NaN) & logMd(165:NaN) & logMGC(1:NaN) & logNgc(6:-inf)check correlated pairs [and interpolate values]
correlation ranking
[x] explore correlations for data type in float
[=] check correlated pairs [and interpolate values]
[=] VMag(418) & KMag(344+74): LR interpolate => KMag(74:NaN), present in errorbar & OLS LR summary
Dist & Ngc(6:0) & Reff(78:NaN)
VMag(418)/KMag(344) & logMd(253)/sigma(271) & logMGC(417)/logNgc(412)
correlation ranking
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.999
Dependent Variable: y AIC: 803.0247
Date: 2018-04-27 12:47 BIC: 807.0602
No. Observations: 418 Log-Likelihood: -400.51
Df Model: 1 F-statistic: 5.005e+05
Df Residuals: 417 Prob (F-statistic): 0.00
R-squared: 0.999 Scale: 0.39886
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 1.1304 0.0016 707.4353 0.0000 1.1273 1.1336
------------------------------------------------------------------
Omnibus: 176.239 Durbin-Watson: 1.879
Prob(Omnibus): 0.000 Jarque-Bera (JB): 1964.668
Skew: 1.479 Prob(JB): 0.000
Kurtosis: 13.201 Condition No.: 1
==================================================================
float[=] check correlated pairs [and interpolate values]
[x] VMag(418) & KMag(344+74)
[=]Dist & Ngc(6:0) & Reff(78:NaN): no significant correlation => leave alone, present in pairplot
VMag(418)/KMag(344) & logMd(253)/sigma(271) & logMGC(417)/logNgc(412)
correlation ranking
float[x] check correlated pairs [and interpolate values]
VMag(418) & KMag(344+74)[x] Dist & Ngc(6:0) & Reff(78:NaN)
[=] VMag(418)/KMag(344) & logMd(253)/sigma(271) & logMGC(417)/logNgc(412)
[=] logNgc(6:-inf) & logMGC(1:NaN): Ngc(6:0) deduct from logMGC (Harris et al. 2011), present in errorbar & OLS LR summary
logNgc(412-1) LR interpolate logMGC(1:NaN)logMGC(418-1-6) LR interpolate Ngc/logNgc(6:0/-inf)VMag/KMag (418/344+74) & logNgc/logMGC (412/417)
logMd(253) & sigma(271) & logNgc/logMGC (412/417)
correlation ranking
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.892
Dependent Variable: y AIC: 954.0297
Date: 2018-04-27 12:47 BIC: 958.0651
No. Observations: 418 Log-Likelihood: -476.01
Df Model: 1 F-statistic: 3437.
Df Residuals: 417 Prob (F-statistic): 1.75e-203
R-squared: 0.892 Scale: 0.57241
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 0.2906 0.0050 58.6292 0.0000 0.2808 0.3003
------------------------------------------------------------------
Omnibus: 40.119 Durbin-Watson: 1.430
Prob(Omnibus): 0.000 Jarque-Bera (JB): 18.367
Skew: 0.316 Prob(JB): 0.000
Kurtosis: 2.191 Condition No.: 1
==================================================================
[x] explore correlations for data type in float
[x] check correlated pairs [and interpolate values]
[x] VMag(418) & KMag(344+74)
[x] Dist & Ngc(412+6) & Reff(78:NaN)
[=] VMag(418)/KMag(344) & logMd(253)/sigma(271) & logMGC(417)/logNgc(412)
[x] logNgc(412+6) & logMGC(417+1)
[=] VMag/KMag (418/344+74) & logNgc/logMGC (412/417): all complete
VMag(original) vs logNgc (more intuitive) logMd(253) & sigma(271) & logNgc/logMGC (412/417)
correlation ranking
[x] explore correlations for data type in float
[x] check correlated pairs [and interpolate values]
[x] VMag(418) & KMag(344+74)
[x] Dist & Ngc(6:0) & Reff(78:NaN)
[=] VMag(418)/KMag(344) & logMd(253)/sigma(271) & logMGC(417+1)/logNgc(412+6)
[x] logNgc(412+6) & logMGC(417+1)
[x] VMag/KMag (418/344+74) & logNgc/logMGC (412+6/417+1)
[=] logMd(253) & sigma(271) & logNgc/logMGC (412+6/417+1)
logMd(253) = f(sigma(271),Reff(340)) vs logNgc(412+6):logMd(253) limiting termcorrelation ranking
[x] explore correlations for data type in float
[x] check correlated pairs [and interpolate values]
[x] VMag & KMag
[x] Dist & Ngc & Reff
[x] VMag(KMag) & logMd (sigma) & logMGC(logNgc)
[=] correlation ranking: Random Forest with feature importance, explain in heatmap
VMag/KMag(344+74) vs logNgc(412+6)/logMGC(417+1)KMag & logMGC > KMag & logNgc > VMag & logNgcKMag & logNgc > KMag & logMGC > VMag & logMGCKMag & logNgc > KMag & logMGC > VMag & logMGCKMag & Ngc > logMd & logMGC > VMag & logMGC, max correlation score dropslogMd vs logMGC/logNgc¶[=] dynamical mass model logMGC/logNgc $=\mathscr{LR}$(logMd)
[=] scatter plot with LR: strong correlation, present in scatter + regplot
OLS LR summary
luminosity model $\mathscr{LR}$(VMag/KMag)
revisit with errorbar
logMd vs logMGC/logNgc¶[=] dynamical mass model logMGC/logNgc $=\mathscr{LR}$(logMd)
[x] scatter plot with LR
[=] OLS LR summary: for logMd(253): logMGC(0.996) > logNgc(0.923)
luminosity model $\mathscr{LR}$(VMag/KMag)
revisit with errorbar
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.923
Dependent Variable: y AIC: 566.4509
Date: 2018-04-27 12:48 BIC: 569.9842
No. Observations: 253 Log-Likelihood: -282.23
Df Model: 1 F-statistic: 3037.
Df Residuals: 252 Prob (F-statistic): 1.41e-142
R-squared: 0.923 Scale: 0.54723
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 0.2379 0.0043 55.1067 0.0000 0.2294 0.2464
------------------------------------------------------------------
Omnibus: 10.559 Durbin-Watson: 1.338
Prob(Omnibus): 0.005 Jarque-Bera (JB): 5.680
Skew: 0.159 Prob(JB): 0.058
Kurtosis: 2.338 Condition No.: 1
==================================================================
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.996
Dependent Variable: y AIC: 340.7403
Date: 2018-04-27 12:48 BIC: 344.2737
No. Observations: 253 Log-Likelihood: -169.37
Df Model: 1 F-statistic: 7.156e+04
Df Residuals: 252 Prob (F-statistic): 2.49e-311
R-squared: 0.996 Scale: 0.22424
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 0.7394 0.0028 267.5077 0.0000 0.7339 0.7448
------------------------------------------------------------------
Omnibus: 5.236 Durbin-Watson: 1.515
Prob(Omnibus): 0.073 Jarque-Bera (JB): 7.433
Skew: -0.017 Prob(JB): 0.024
Kurtosis: 3.839 Condition No.: 1
==================================================================
VMag/KMag vs logNgc¶[x] dynamical mass model $\mathscr{LR}$(logMd)
[=] luminosity model logNgc $=\mathscr{LR}$(VMag/KMag)
[=] scatter plot with LR: strong correlation, present in scatter + regplot
OLS LR summary
revisit with errorbar
VMag/KMag vs logNgc¶[x] dynamical mass model $\mathscr{LR}$(logMd)
[x] luminosity model logNgc=$\mathscr{LR}$(VMag/KMag)
[x] scatter plot with LR
[=] OLS LR summary
VMag(0.866) < KMag(0.869) < logMd(0.923/0.996)logMd(253/418), KMag(344/418) interpolated from VMaglogMd requires good photometry dataKMagrevisit
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.866
Dependent Variable: y AIC: 1043.0025
Date: 2018-04-27 12:48 BIC: 1047.0379
No. Observations: 418 Log-Likelihood: -520.50
Df Model: 1 F-statistic: 2698.
Df Residuals: 417 Prob (F-statistic): 3.34e-184
R-squared: 0.866 Scale: 0.70819
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 -0.1106 0.0021 -51.9461 0.0000 -0.1148 -0.1064
------------------------------------------------------------------
Omnibus: 24.883 Durbin-Watson: 1.467
Prob(Omnibus): 0.000 Jarque-Bera (JB): 19.434
Skew: 0.434 Prob(JB): 0.000
Kurtosis: 2.397 Condition No.: 1
==================================================================
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.870
Dependent Variable: y AIC: 1030.4039
Date: 2018-04-27 12:48 BIC: 1034.4393
No. Observations: 418 Log-Likelihood: -514.20
Df Model: 1 F-statistic: 2794.
Df Residuals: 417 Prob (F-statistic): 6.21e-187
R-squared: 0.870 Scale: 0.68716
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 -0.0980 0.0019 -52.8557 0.0000 -0.1017 -0.0944
------------------------------------------------------------------
Omnibus: 22.487 Durbin-Watson: 1.448
Prob(Omnibus): 0.000 Jarque-Bera (JB): 19.376
Skew: 0.453 Prob(JB): 0.000
Kurtosis: 2.458 Condition No.: 1
==================================================================
[x] dynamical mass model logMGC/logNgc=$\mathscr{LR}$(logMd)
[x] luminosity model logNgc=$\mathscr{LR}$(VMag/KMag)
[=] revisit with errorbar, present in 2 errorbars
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.870
Dependent Variable: y AIC: 1030.4039
Date: 2018-04-27 12:48 BIC: 1034.4393
No. Observations: 418 Log-Likelihood: -514.20
Df Model: 1 F-statistic: 2794.
Df Residuals: 417 Prob (F-statistic): 6.21e-187
R-squared: 0.870 Scale: 0.68716
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 -0.0980 0.0019 -52.8557 0.0000 -0.1017 -0.0944
------------------------------------------------------------------
Omnibus: 22.487 Durbin-Watson: 1.448
Prob(Omnibus): 0.000 Jarque-Bera (JB): 19.376
Skew: 0.453 Prob(JB): 0.000
Kurtosis: 2.458 Condition No.: 1
==================================================================
logMd vs logMGC/logNgc and VMag/KMag vs logNgc¶[x] dynamical mass model logMGC/logNgc=$\mathscr{LR}$(logMd)
[x] luminosity model logNgc=$\mathscr{LR}$(VMag/KMag)
[=] revisit with errorbar, present in 2 errorbars
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.870
Dependent Variable: y AIC: 1030.4039
Date: 2018-04-27 12:48 BIC: 1034.4393
No. Observations: 418 Log-Likelihood: -514.20
Df Model: 1 F-statistic: 2794.
Df Residuals: 417 Prob (F-statistic): 6.21e-187
R-squared: 0.870 Scale: 0.68716
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 -0.0980 0.0019 -52.8557 0.0000 -0.1017 -0.0944
------------------------------------------------------------------
Omnibus: 22.487 Durbin-Watson: 1.448
Prob(Omnibus): 0.000 Jarque-Bera (JB): 19.376
Skew: 0.453 Prob(JB): 0.000
Kurtosis: 2.458 Condition No.: 1
==================================================================
VMag, normalized to $M_V$ = −15, present in scatter[=] LOWESS regression for $S_N$ vs VMag, present in plot
error comparison
Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.870
Dependent Variable: y AIC: 1030.4039
Date: 2018-04-27 12:48 BIC: 1034.4393
No. Observations: 418 Log-Likelihood: -514.20
Df Model: 1 F-statistic: 2794.
Df Residuals: 417 Prob (F-statistic): 6.21e-187
R-squared: 0.870 Scale: 0.68716
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 -0.0980 0.0019 -52.8557 0.0000 -0.1017 -0.0944
------------------------------------------------------------------
Omnibus: 22.487 Durbin-Watson: 1.448
Prob(Omnibus): 0.000 Jarque-Bera (JB): 19.376
Skew: 0.453 Prob(JB): 0.000
Kurtosis: 2.458 Condition No.: 1
==================================================================
VMag[x] LOWESS regression for $S_N$ vs VMag
[=] error comparison, zoom in view, justificationin errorbar
VMag range [-11.17, -24.19]: lower-bound Results: Ordinary least squares
==================================================================
Model: OLS Adj. R-squared: 0.870
Dependent Variable: y AIC: 1030.4039
Date: 2018-04-27 12:48 BIC: 1034.4393
No. Observations: 418 Log-Likelihood: -514.20
Df Model: 1 F-statistic: 2794.
Df Residuals: 417 Prob (F-statistic): 6.21e-187
R-squared: 0.870 Scale: 0.68716
---------------------------------------------------------------------
Coef. Std.Err. t P>|t| [0.025 0.975]
---------------------------------------------------------------------
x1 -0.0980 0.0019 -52.8557 0.0000 -0.1017 -0.0944
------------------------------------------------------------------
Omnibus: 22.487 Durbin-Watson: 1.448
Prob(Omnibus): 0.000 Jarque-Bera (JB): 19.376
Skew: 0.453 Prob(JB): 0.000
Kurtosis: 2.458 Condition No.: 1
==================================================================
[=] bandpass conversion: B/I/K in GWGC from White et al. (2011) to V/K in GCSN model from Harris et al. (2013)
[=] GWGC luminosities, Abs_Mag_B/I(7,877/5,606), App_Mag_B/I/K(7,877/5,606/2,995), and err_Abs/App_Mag_B/I/K:
explore correlations
fit conversion constant
choose most probable value
revert to Ngc
GC distribution
| App_Mag_B | err_App_Mag_B | Abs_Mag_B | err_Abs_Mag_B | App_Mag_I | err_App_Mag_I | Abs_Mag_I | err_Abs_Mag_I | App_Mag_K | err_App_Mag_K | |
|---|---|---|---|---|---|---|---|---|---|---|
| count | 7877.000000 | 7877.000000 | 7877.000000 | 7877.000000 | 5606.000000 | 5606.000000 | 5606.000000 | 0.0 | 2995.000000 | 2287.00000 |
| min | -4.600000 | 0.040000 | -22.590000 | 0.070000 | 5.730000 | 0.030000 | -25.520000 | NaN | 7.271000 | 0.02000 |
| max | 25.260000 | 3.020000 | -4.970000 | 3.020000 | 30.440000 | 6.880000 | -1.680000 | NaN | 15.978000 | 0.29000 |
| mean | 14.981126 | 0.377188 | -16.210734 | 0.387227 | 14.158673 | 0.174930 | -17.158776 | NaN | 12.968526 | 0.10208 |
| std | 2.280469 | 0.176294 | 2.451796 | 0.175116 | 2.463852 | 0.317297 | 2.627928 | NaN | 1.594678 | 0.04896 |
| NaNs | 1069.000000 | 1069.000000 | 1069.000000 | 1069.000000 | 3340.000000 | 3340.000000 | 3340.000000 | 8946.0 | 5951.000000 | 6659.00000 |
[=] bandpass conversion: B/I/K in GWGC from White et al. (2011) to V/K in GCSN model from Harris et al. (2013)
[x] GWGC luminosities
[=] explore correlations, explain in pairplots with heatmap
xxx & err_xxx: intuitive, not interestingabs_x & app_x: b'cos $app\_x-abs\_X=5\log_{10}(d)-5$err_Abs_x & err_App_x: of course, passfit conversion constant
choose most probable value
revert to Ngc
GC distribution
[=] bandpass conversion: B/I/K in GWGC from White et al. (2011) to V/K in GCSN model from Harris et al. (2013)
[x] GWGC luminosities
[x] explore correlations
[=] fit conversion constant, 197 galaxies in both catalogs, inspect in mag/color scatter
[Abs_Mag_B (7,877), err_Abs_Mag_B (7,877)] > [Abs_Mag_I (5,606), err_Abs_Mag_I (0)]: => prefer B over I for Ngc estimation B/I in GWGC (197/148) vs V/K in GCSN=> almost homogenous, prefer B for fitB/I in GWGC (197/148) vs V-K in GCSN=> color/Mcat didn't help (unequal axis)choose most probable value
revert to Ngc
GC distribution
[=] bandpass conversion: B/I/K in GWGC from White et al. (2011) to V/K in GCSN model from Harris et al. (2013)
[x] GWGC luminosities
[x] explore correlations
[x] fit conversion constant: B/I in GWGC (197/148) vs V/K in GC SN on 197 shared galaxies
[=] choose most probable value, based on kdeplot
revert to Ngc
GC distribution
[x] bandpass conversion: B/I/K in GWGC from White et al. (2011) to V/K in GCSN model from Harris et al. (2013)
[x] GWGC luminosities
[x] explore correlations
[x] fit conversion constant: B/I in GWGC (197/148) vs V/K in GC SN on 197 shared galaxies
[x] choose most probable value, based on kdeplot
[=] revert to Ngc:
BMag/IMag in GWGC to VMag in GCSN: B to V (7,877) + I to V (extra 32) = 7,909/8,946658,910, over 8,946 galaxies within 30 MpcmodestGC distribution
| Name | RA | Dec | VMag | Dist | Abs_Mag_B | Abs_Mag_I | Ngc | e_Ngc | iMType | |
|---|---|---|---|---|---|---|---|---|---|---|
| count | 8946 | 8946.000000 | 8946.000000 | 7909.000000 | 8946.000000 | 7877.000000 | 5606.000000 | 8946.000000 | 8946.000000 | 8946.000000 |
| min | NaN | 0.005940 | -89.334590 | -23.478214 | 0.010000 | -22.590000 | -25.520000 | 0.000000 | 0.000000 | 0.000000 |
| max | NaN | 23.996680 | 86.131800 | -5.858214 | 29.986000 | -4.970000 | -1.680000 | 11608.000000 | 7187.308776 | 6.000000 |
| mean | NaN | 10.973250 | 3.455383 | -17.094336 | 19.011495 | -16.210734 | -17.158776 | 74.085848 | 46.053841 | 0.916387 |
| std | NaN | 5.277991 | 34.889882 | 2.452026 | 6.616150 | 2.451796 | 2.627928 | 262.963677 | 190.324592 | 0.936047 |
| 168135 | SDSSJ142852.74+123454.0 | 14.481330 | 12.581660 | NaN | 28.958000 | NaN | NaN | 0.000000 | 0.000000 | 0.000000 |
| 168242 | SDSSJ144213.99+333517.4 | 14.703880 | 33.588320 | NaN | 24.528000 | NaN | NaN | 0.000000 | 0.000000 | 0.000000 |
| 146627 | SDSSJ080937.39+413526.6 | 8.160380 | 41.590710 | -14.168214 | 11.625000 | -13.280000 | -12.870000 | 5.000000 | 2.037453 | 0.000000 |
| 24031 | UGC06251 | 11.223940 | 53.595040 | -17.718214 | 18.450000 | -16.830000 | NaN | 25.000000 | 7.388049 | 1.000000 |
| 171592 | SDSSJ151605.06+561614.7 | 15.268080 | 56.270820 | -16.628214 | 9.875000 | -15.740000 | -16.790000 | 14.000000 | 3.364954 | 0.000000 |
[x] bandpass conversion: B/I/K in GWGC from White et al. (2011) to V/K in GCSN model from Harris et al. (2013)
[x] GWGC luminosities
[x] explore correlations
[x] fit conversion constant: B/I in GWGC (197/148) vs V/K in GC SN on 197 shared galaxies
[x] choose most probable value, based on kdeplot
[x] revert to Ngc: 7,909/8,946 with Ngc
658,910, over 8,946 galaxies within 30 Mpcmodest[=] GC distribution
| Name | RA | Dec | VMag | Dist | Abs_Mag_B | Abs_Mag_I | Ngc | e_Ngc | iMType | |
|---|---|---|---|---|---|---|---|---|---|---|
| count | 8946 | 8946.000000 | 8946.000000 | 0.0 | 8946.000000 | 7877.000000 | 5606.000000 | 0.0 | 0.0 | 8946.000000 |
| min | NaN | 0.005940 | -89.334590 | NaN | 0.010000 | -22.590000 | -25.520000 | NaN | NaN | 0.000000 |
| max | NaN | 23.996680 | 86.131800 | NaN | 29.986000 | -4.970000 | -1.680000 | NaN | NaN | 6.000000 |
| mean | NaN | 10.973250 | 3.455383 | NaN | 19.011495 | -16.210734 | -17.158776 | NaN | NaN | 0.916387 |
| std | NaN | 5.277991 | 34.889882 | NaN | 6.616150 | 2.451796 | 2.627928 | NaN | NaN | 0.936047 |
| 168135 | SDSSJ142852.74+123454.0 | 14.481330 | 12.581660 | NaN | 28.958000 | NaN | NaN | NaN | NaN | 0.000000 |
| 168242 | SDSSJ144213.99+333517.4 | 14.703880 | 33.588320 | NaN | 24.528000 | NaN | NaN | NaN | NaN | 0.000000 |
| 146627 | SDSSJ080937.39+413526.6 | 8.160380 | 41.590710 | NaN | 11.625000 | -13.280000 | -12.870000 | NaN | NaN | 0.000000 |
| 24031 | UGC06251 | 11.223940 | 53.595040 | NaN | 18.450000 | -16.830000 | NaN | NaN | NaN | 1.000000 |
| 171592 | SDSSJ151605.06+561614.7 | 15.268080 | 56.270820 | NaN | 9.875000 | -15.740000 | -16.790000 | NaN | NaN | 0.000000 |
General setup
Variations of parameters
Age spread: the Hubble time is assumed to be 13.5 Gyr
GC simulations on TACC
GClib LIBRARY: TO BE RELEASED, 652,196 -> 662,772


Assembling BBH database
BBHlib: 87,365 over 3,240 GCs
LoGal: 7,909 galaxies with GCs within 30 Mpc
GClib: 662,772 GCs over 7,909 galaxies
=> BBHdata: 17,883,760 BBHs


Assembling BBH database
BBHlib: extracted BBHs from 3,240 GC simulations
LoGal: galaxy within 30 Mpc
GClib: 662,772 GCs over 7,909 Galaxy
assign random GC model
GClib model diversity:
=> BBHdata: 17,883,760 BBHs
```
for each GC in LoGal.Ngc:
assign random GC model
get galaxy info. from LoGal
get BBH from BBHlib with GC model
combine
```
BBHdata
| host galaxy info | GC info | BBH info |
|---|---|---|
| 1-5 | 6-7 | 8-14 |
| GC, RA, DEC, Dist, VMag | Model, T_GC (13.5Gyr-Age)[Fe/H, f_b, ...] | T_eject, M1, M2, Seperation,Ecc, (BBHlib.Model=GClib.Model) |
| T_eject | M1 | M2 | Seperation | Ecc | Model | |
|---|---|---|---|---|---|---|
| 0 | 0.0079 | 26.460 | 11.6190 | 36.1060 | 0.3173 | 187-7 |
| 1 | 0.0118 | 22.999 | 7.8955 | 222.4900 | 0.0000 | 187-7 |
| 2 | 0.0157 | 23.190 | 8.1896 | 20.0940 | 0.0978 | 187-7 |
| 3 | 0.0315 | 23.175 | 8.2358 | 275.1300 | 0.4285 | 187-7 |
| 4 | 0.0315 | 22.895 | 7.7428 | 91.2530 | 0.8577 | 187-7 |
| 5 | 0.0315 | 23.874 | 9.4294 | 173.7900 | 0.7756 | 187-7 |
| 6 | 0.0315 | 25.761 | 10.9290 | 868.3400 | 0.0000 | 187-7 |
| 7 | 0.0499 | 25.777 | 11.1890 | 793.2600 | 0.3369 | 187-7 |
| 8 | 0.4144 | 13.386 | 25.7790 | 78.5730 | 0.8045 | 187-7 |
| 9 | 0.7613 | 26.697 | 26.9370 | 185.3200 | 0.5458 | 187-7 |
| 10 | 0.7886 | 37.914 | 28.3050 | 331.4400 | 0.3656 | 187-7 |
| 11 | 0.8329 | 26.815 | 26.1680 | 191.4900 | 0.8724 | 187-7 |
| 12 | 1.0465 | 23.922 | 24.4750 | 51.2410 | 0.6847 | 187-7 |
| 13 | 1.2962 | 25.195 | 25.3110 | 217.5200 | 0.5197 | 187-7 |
| 14 | 1.5873 | 26.058 | 23.4010 | 299.2500 | 0.9600 | 187-7 |
| 15 | 1.6343 | 26.829 | 26.8670 | 738.9600 | 0.7871 | 187-7 |
| 16 | 2.1382 | 25.227 | 26.1970 | 203.5100 | 0.6862 | 187-7 |
| 17 | 2.2051 | 26.931 | 24.2770 | 94.9120 | 0.4887 | 187-7 |
| 18 | 2.2264 | 26.946 | 24.3750 | 334.5500 | 0.7257 | 187-7 |
| 19 | 2.2484 | 27.270 | 26.6270 | 317.5200 | 0.7703 | 187-7 |
| 20 | 2.6205 | 23.248 | 25.5000 | 567.8200 | 0.9503 | 187-7 |
| 21 | 2.9080 | 21.622 | 24.2440 | 306.3300 | 0.5750 | 187-7 |
| 22 | 2.9240 | 23.420 | 23.5020 | 623.8000 | 0.3999 | 187-7 |
| 23 | 2.9433 | 23.193 | 40.5950 | 249.4500 | 0.8743 | 187-7 |
| 24 | 3.0288 | 23.418 | 22.5560 | 321.2600 | 0.8533 | 187-7 |
| 25 | 3.1998 | 26.911 | 23.5720 | 706.2900 | 0.6283 | 187-7 |
| 26 | 3.3023 | 22.814 | 24.3720 | 117.9900 | 0.9995 | 187-7 |
| 27 | 3.5353 | 22.164 | 20.8490 | 175.1500 | 0.9920 | 187-7 |
| 28 | 4.0062 | 13.906 | 13.7430 | 9.3148 | 0.8166 | 187-7 |
| 29 | 4.1225 | 25.177 | 26.7480 | 131.3900 | 0.2675 | 187-7 |
| 30 | 4.1471 | 25.910 | 52.3300 | 326.6700 | 0.6479 | 187-7 |
| 31 | 4.2271 | 24.400 | 13.2930 | 12.6890 | 0.6893 | 187-7 |
| 32 | 4.2399 | 23.670 | 24.8380 | 566.1000 | 0.8899 | 187-7 |
| 33 | 4.3531 | 20.876 | 28.1520 | 428.8500 | 0.5484 | 187-7 |
| 34 | 5.4659 | 24.859 | 22.3790 | 292.1900 | 0.9335 | 187-7 |
| 35 | 5.4829 | 23.489 | 21.2020 | 133.5700 | 0.6346 | 187-7 |
| 37 | 5.8268 | 24.884 | 49.2470 | 383.3000 | 0.6541 | 187-7 |
| 38 | 6.0204 | 27.398 | 19.7820 | 188.6200 | 0.4226 | 187-7 |
| 39 | 6.3426 | 22.117 | 23.3930 | 79.8000 | 0.4630 | 187-7 |
| 40 | 6.6529 | 27.150 | 11.6090 | 81.0850 | 0.7587 | 187-7 |
| 41 | 7.0097 | 22.594 | 27.0070 | 120.4400 | 0.4160 | 187-7 |
| 42 | 7.8008 | 24.080 | 26.8740 | 114.6100 | 0.8300 | 187-7 |
| 43 | 7.8008 | 24.624 | 25.3410 | 16084.0000 | 0.8807 | 187-7 |
| 44 | 8.3006 | 13.503 | 18.4360 | 148.1200 | 0.6899 | 187-7 |
| 45 | 8.3461 | 20.889 | 26.7890 | 189.2700 | 0.4582 | 187-7 |
| 47 | 9.4491 | 13.858 | 27.8840 | 138.2300 | 0.6108 | 187-7 |
| 48 | 9.7201 | 22.506 | 22.7170 | 197.7200 | 0.8736 | 187-7 |
| 49 | 10.4762 | 13.017 | 12.7020 | 175.4100 | 0.8080 | 187-7 |
| 50 | 10.7644 | 27.078 | 23.0610 | 243.2300 | 0.5025 | 187-7 |
| 51 | 11.0249 | 11.821 | 13.9440 | 303.7100 | 0.7290 | 187-7 |
| 53 | 11.7761 | 13.669 | 11.8270 | 233.4800 | 0.8816 | 187-7 |
| T_eject | M1 | M2 | Seperation | Ecc | |
|---|---|---|---|---|---|
| count | 87364.000000 | 87364.000000 | 87364.000000 | 8.736400e+04 | 87364.000000 |
| mean | 1.822836 | 24.921859 | 24.182866 | 5.255124e+02 | 0.685132 |
| std | 2.452034 | 17.785691 | 19.599724 | 2.390284e+04 | 0.279524 |
| min | 0.005200 | 3.000100 | 3.014100 | 0.000000e+00 | 0.000000 |
| 25% | 0.237200 | 19.756000 | 13.493000 | 4.714150e+01 | 0.518700 |
| 50% | 0.818000 | 23.190000 | 22.680000 | 9.762050e+01 | 0.765700 |
| 75% | 2.389600 | 26.121000 | 27.013500 | 2.093825e+02 | 0.917400 |
| max | 13.095300 | 2669.500000 | 820.870000 | 5.894200e+06 | 1.000000 |
| Galaxy | RA | Dec | Dist | VMag | Model | Age | |
|---|---|---|---|---|---|---|---|
| 133 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 |
| Galaxy | RA | Dec | Dist | VMag | Age | |
|---|---|---|---|---|---|---|
| count | 662772.000000 | 662772.000000 | 662772.000000 | 662772.000000 | 662772.000000 | 662772.000000 |
| mean | 30603.041664 | 10.959227 | -0.355427 | 21.170577 | -20.534914 | 11.890659 |
| std | 27284.737734 | 5.445347 | 38.155540 | 5.926124 | 1.693741 | 0.864954 |
| min | 25.000000 | 0.007920 | -89.334590 | 0.010000 | -23.478214 | 8.000000 |
| 25% | 15238.000000 | 7.451130 | -28.805970 | 17.783000 | -21.698214 | 11.307087 |
| 50% | 27130.000000 | 11.850490 | -0.480240 | 21.478000 | -20.888214 | 11.968504 |
| 75% | 35109.000000 | 13.284170 | 25.500910 | 25.351000 | -19.788214 | 12.535433 |
| max | 183326.000000 | 23.996680 | 86.131800 | 29.986000 | -11.178214 | 13.480315 |
| Galaxy | RA | Dec | Dist | VMag | Model | Age | T_eject | M1 | M2 | Seperation | Ecc | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0079 | 26.460 | 11.6190 | 36.1060 | 0.3173 |
| 1 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0118 | 22.999 | 7.8955 | 222.4900 | 0.0000 |
| 2 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0157 | 23.190 | 8.1896 | 20.0940 | 0.0978 |
| 3 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0315 | 23.175 | 8.2358 | 275.1300 | 0.4285 |
| 4 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0315 | 22.895 | 7.7428 | 91.2530 | 0.8577 |
| 5 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0315 | 23.874 | 9.4294 | 173.7900 | 0.7756 |
| 6 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0315 | 25.761 | 10.9290 | 868.3400 | 0.0000 |
| 7 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.0499 | 25.777 | 11.1890 | 793.2600 | 0.3369 |
| 8 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.4144 | 13.386 | 25.7790 | 78.5730 | 0.8045 |
| 9 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.7613 | 26.697 | 26.9370 | 185.3200 | 0.5458 |
| 10 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.7886 | 37.914 | 28.3050 | 331.4400 | 0.3656 |
| 11 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 0.8329 | 26.815 | 26.1680 | 191.4900 | 0.8724 |
| 12 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 1.0465 | 23.922 | 24.4750 | 51.2410 | 0.6847 |
| 13 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 1.2962 | 25.195 | 25.3110 | 217.5200 | 0.5197 |
| 14 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 1.5873 | 26.058 | 23.4010 | 299.2500 | 0.9600 |
| 15 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 1.6343 | 26.829 | 26.8670 | 738.9600 | 0.7871 |
| 16 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.1382 | 25.227 | 26.1970 | 203.5100 | 0.6862 |
| 17 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.2051 | 26.931 | 24.2770 | 94.9120 | 0.4887 |
| 18 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.2264 | 26.946 | 24.3750 | 334.5500 | 0.7257 |
| 19 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.2484 | 27.270 | 26.6270 | 317.5200 | 0.7703 |
| 20 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.6205 | 23.248 | 25.5000 | 567.8200 | 0.9503 |
| 21 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.9080 | 21.622 | 24.2440 | 306.3300 | 0.5750 |
| 22 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.9240 | 23.420 | 23.5020 | 623.8000 | 0.3999 |
| 23 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 2.9433 | 23.193 | 40.5950 | 249.4500 | 0.8743 |
| 24 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 3.0288 | 23.418 | 22.5560 | 321.2600 | 0.8533 |
| 25 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 3.1998 | 26.911 | 23.5720 | 706.2900 | 0.6283 |
| 26 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 3.3023 | 22.814 | 24.3720 | 117.9900 | 0.9995 |
| 27 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 3.5353 | 22.164 | 20.8490 | 175.1500 | 0.9920 |
| 28 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 4.0062 | 13.906 | 13.7430 | 9.3148 | 0.8166 |
| 29 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 4.1225 | 25.177 | 26.7480 | 131.3900 | 0.2675 |
| 30 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 4.1471 | 25.910 | 52.3300 | 326.6700 | 0.6479 |
| 31 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 4.2271 | 24.400 | 13.2930 | 12.6890 | 0.6893 |
| 32 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 4.2399 | 23.670 | 24.8380 | 566.1000 | 0.8899 |
| 33 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 4.3531 | 20.876 | 28.1520 | 428.8500 | 0.5484 |
| 34 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 5.4659 | 24.859 | 22.3790 | 292.1900 | 0.9335 |
| 35 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 5.4829 | 23.489 | 21.2020 | 133.5700 | 0.6346 |
| 36 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 5.8268 | 24.884 | 49.2470 | 383.3000 | 0.6541 |
| 37 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 6.0204 | 27.398 | 19.7820 | 188.6200 | 0.4226 |
| 38 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 6.3426 | 22.117 | 23.3930 | 79.8000 | 0.4630 |
| 39 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 6.6529 | 27.150 | 11.6090 | 81.0850 | 0.7587 |
| 40 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 7.0097 | 22.594 | 27.0070 | 120.4400 | 0.4160 |
| 41 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 7.8008 | 24.080 | 26.8740 | 114.6100 | 0.8300 |
| 42 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 7.8008 | 24.624 | 25.3410 | 16084.0000 | 0.8807 |
| 43 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 8.3006 | 13.503 | 18.4360 | 148.1200 | 0.6899 |
| 44 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 8.3461 | 20.889 | 26.7890 | 189.2700 | 0.4582 |
| 45 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 9.4491 | 13.858 | 27.8840 | 138.2300 | 0.6108 |
| 46 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 9.7201 | 22.506 | 22.7170 | 197.7200 | 0.8736 |
| 47 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 10.4762 | 13.017 | 12.7020 | 175.4100 | 0.8080 |
| 48 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 10.7644 | 27.078 | 23.0610 | 243.2300 | 0.5025 |
| 49 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 11.0249 | 11.821 | 13.9440 | 303.7100 | 0.7290 |
| 50 | 133 | 0.05414 | 16.14555 | 13.183 | -20.568214 | 187-7 | 11.070866 | 11.7761 | 13.669 | 11.8270 | 233.4800 | 0.8816 |
| PGC | RA | Dec | Dist | VMag | Age | T_eject | M1 | M2 | Seperation | Ecc | Period | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| count | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 1.788376e+07 | 0.0 |
| mean | 3.059303e+04 | 1.095810e+01 | -4.172814e-01 | 2.116666e+01 | -2.053641e+01 | 1.189068e+01 | 1.823276e+00 | 2.492480e+01 | 2.418357e+01 | 5.184244e+02 | 6.851611e-01 | NaN |
| std | 2.726805e+04 | 5.447050e+00 | 3.812537e+01 | 5.926283e+00 | 1.693210e+00 | 8.648354e-01 | 2.451835e+00 | 1.776007e+01 | 1.958754e+01 | 2.311852e+04 | 2.795283e-01 | NaN |
| min | 2.500000e+01 | 7.920000e-03 | -8.933459e+01 | 1.000000e-02 | -2.347821e+01 | 8.000000e+00 | 5.200000e-03 | 3.000100e+00 | 3.014100e+00 | 0.000000e+00 | 0.000000e+00 | NaN |
| 25% | 1.523800e+04 | 7.451130e+00 | -2.880597e+01 | 1.778300e+01 | -2.169821e+01 | 1.130709e+01 | 2.375000e-01 | 1.975600e+01 | 1.349300e+01 | 4.715300e+01 | 5.188000e-01 | NaN |
| 50% | 2.713300e+04 | 1.185049e+01 | -4.986800e-01 | 2.147800e+01 | -2.088821e+01 | 1.196850e+01 | 8.189000e-01 | 2.319300e+01 | 2.268100e+01 | 9.761800e+01 | 7.658000e-01 | NaN |
| 75% | 3.510900e+04 | 1.328468e+01 | 2.550091e+01 | 2.535100e+01 | -1.978821e+01 | 1.253543e+01 | 2.389600e+00 | 2.612100e+01 | 2.700700e+01 | 2.093100e+02 | 9.175000e-01 | NaN |
| max | 1.833260e+05 | 2.399668e+01 | 8.613180e+01 | 2.998600e+01 | -1.117821e+01 | 1.348031e+01 | 1.309530e+01 | 2.669500e+03 | 8.208700e+02 | 5.894200e+06 | 1.000000e+00 | NaN |
formulae: Peters (1964)
code: Parallelized from Benacquista (2004)
Separation, Ecc, PeriodMerger event rate (6.1.3)




contact: contact@domij.info
github repository: LocalBBH
Resume: 03/24/2018 cache